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Abstract: 

We discuss corrections to the ratio of shear viscosity to entropy density rj/ s in higher- 
derivative gravity theories. Generically, these theories contain ghost modes with 
Planck-scale masses. Motivated by general considerations about unitarity, we pro- 
pose new boundary conditions for the equations of motion of the graviton perturba- 
tions that force the amplitude of the ghosts modes to vanish. We analyze explicitly 
four-derivative perturbative corrections to Einstein gravity which generically lead to 
four-derivative equations of motion, compare our choice of boundary conditions to 
previous proposals and show that, with our new prescription, the ratio rj/ s remains 
at the Einstein-gravity value of l/47r to leading order in the corrections. It is argued 
that, when the new boundary conditions are imposed on six and higher-derivative 
equations of motion, rj/ s can only increase from the Einstein-gravity value. We also 
recall some general arguments that support the validity of our results to all orders 
in the strength of the corrections to Einstein gravity. We then discuss the partic- 
ular case of Gauss-Bonnet gravity, for which the equations of motion are only of 
two-derivative order and the value of r]/ s can decrease below l/47r when treated in 
a nonperturbative way. Our findings provide further evidence for the validity of the 
KSS bound for theories that can be viewed as perturbative corrections to Einstein 
Gravity. 
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1. Introduction 

The gauge-gravity duality [|l], |^ allows one to holographically map black brane ther- 
modynamics and hydrodynamics in the Anti-deSitter (AdS) bulk to their gauge- 
theory correspondents at the AdS boundary O, M . 
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For a large class of strongly coupled fluids (essentially, any with a two-derivative 
gravity dual), the ratio of the shear viscosity t] to the entropy density s is, in ap- 



propriately chosen units, remarkably low: f = ^ H) ^- The uncertainty principle 
can be used to argue that rj/s should have a lower bound of order unity P] ; leading 
Kovtun, Son and Starinets (KSS) to propose that I/Att is a universal lower bound on 
this ratio. However, investigations of higher- derivative gravity theories have revealed 
that the so-called KSS bound can apparently be violated [0, §]. Many subsequent in- 
quiries have appealed to gauge-theory causality as the physical principle that bounds 
v/s (e.g.) §-0. 

We have recently asserted |TH] that a unitary ghost-free extension of Einstein 



gravity can, at best, saturate the KSS bound. We have observed that 77 is a grav- 
itational coupling and any such coupling can only increase from its Einstein value 
in unitary theories |TB|, |T^ |IB|. The source of tension between our claim and the 



apparently bound- violating gravity theories is that, even though higher-derivative 
gravity theories contain ghosts, these are typically at the Planck scale and therefore 
expected to be irrelevant to the calculation of hydrodynamic transport coefficients 
like f]. 

Let us briefly recall why ghosts are problematic for any classical or quantum 
field theory (also see, e.g., [|l^]-p3|). If ghosts are present, then any state, including 
the vacuum, would be catastrophically unstable due to the spontaneous creation of 
positive and negative energy particles having zero total energy. For a classical theory 
of gravity, in particular, ghost gravitons lead to instabilities at arbitrarily short time 
scales or superluminal propagation in both the gravitational and matter sectors. 
These issues have to do with the wrong sign of the kinetic energy and not with the 
relative sign of the kinetic and mass terms. Although the latter can, if "incorrectly" 
chosen, have further effects that undermine both stability and causality. 

The necessity to remove ghosts becomes even more acute when the gauge-gravity 
duality is considered. From the field-theory side, only a finite number of states is 
permissible, which corresponds to a truncated spectrum of bulk excitations. This 
so-called "stringy exclusion principle" [E3 has been argued to prohibit ghosts 
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From an effective-theory perspective, the situation is not so clear cut when all 
the ghosts have a mass on the order of the ultraviolet cutoff. One might then reason 
that it is now safe to disregard their presence. This could be true in many instances, 
but this reasoning will be shown to be incorrect for a calculation of rj. So, to proceed 
in a sensible manner, one should ensure that the ghost has decoupled before the 
computation is carried out. This could be accomplished through a choice of boundary 
conditions (BC's). 

Later on, we use four- derivative theories of gravity to explicitly show that the 
Planck-scale ghosts do indeed infiltrate the existing schemes of calculating rj. Then, 
motivated by considerations of unitarity, we propose a new set of BC's for the higher- 
order equations of motion that forces the amplitude of the ghost modes to vanish. 
Applying the new prescription to four-derivative gravity theories, we find that rj/s = 
1/Att to leading order in the strength of the corrections to Einstein gravity. We then 
argue that, after imposing the new BC's on six and higher-derivative equations of 
motion, 7]/s can only increase above 1/4it. These results support our previous claim 
that l/47r is indeed a lower bound on rj/s for any unitary weakly coupled extension 
of Einstein gravity. 

Let us now suppose that a gauge field theory is given and its AdS gravity dual is 
known. So, just how does one go about determining rj/s from a black brane theory 
on the AdS side? For any two-derivative theory of gravity, ^ this is a well understood 
matter, and the resulting answer of any method must agree with the well-known 
Einstein result of rj/s = l/Air [^, Q. But, for higher-derivative gravity theories, the 
situation is still quite ambiguous. 

The Wald Noether-charge formalism [^U], |2^ is appropriate and accepted for the 



calculation of the entropy or its density s. For the shear viscosity rj, the situation is 
less clear. The earliest calculations of rj for Einstein gravity used holographic tech- 
niques to express the field-theory viscosity — which follows from the Kubo formula 
— in terms of horizon- valued parameters of the AdS brane theory. (For a review and 



^The two-derivative class is meant to include theories that are effectively two-derivative; e.g.. 
any f{TZ) theory can be expressed as Einstein gravity plus a scalar. 
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references, see [§.) Later, the physical connection between the boundary and black 
brane theories has been made more explicit via the membrane-fluid interpretation of 
horizon hydrodjTiamics 0. 

For any two-derivative gravity theory, one can also obtain rj from the horizon 
value of the coefficient of the kinetic term of the transversely polarized gravitons 



[Q-[^. This is the same as reading off the coefficient of the graviton propagator 
(hxyhxy), where h^^, is the first-order correction to the background metric ^. 

For higher- derivative theories, it is reasonable to apply a similar procedure. Ex- 
cept that, in this case, one first requires a well-defined prescription for dealing with 
the higher (than second) order derivatives which inevitably turn up in the action. 
A popular approach has appeared in different yet equivalent guises in the literature 
[0-[0. For instance, according to [Q, one is instructed to iteratively apply the 



field equation to the action, consistently trading off the perturbative-order higher 
derivatives for lower ones (while ignoring terms of second or higher perturbative 
order), until a two- derivative action has been obtained. Then, the viscosity is ex- 
tracted from the kinetic term of this effective action. Identical results have been 
obtained from a canonical-momentum formulation of rj [^ , which can be viewed as 



a higher-derivative generalization of the membrane paradigm 0] as interpreted in 

An alternative approach [^ prescribes ignoring the higher-order derivatives and 
extracting t] straight from the kinetic term. In this approach, 77 is presumed to 
maintain its identity as the coefficient of the propagator. For theories limited to four 
derivatives, the two approaches are in agreement. This agreement does not, however, 
persist for theories with six or more derivatives, nor could it be expected to. For 
the sake of completeness, we will elaborate further on these methods in an appendix. 
Our goal in this paper is not to argue for the merit of one prescription versus the 
other. Rather, it is our contention that both of these approaches should be corrected. 

The rest of the paper proceeds as follows: In the next section, we introduce 



^x and y represent transverse brane directions that are mutually orthogonal to one another, as 
well as to the propagating direction of gravitons moving along the brane. 



-4- 



the explicit models which are four- derivative extensions of Einstein gravity. Here, 
our conventions are set and much of the necessary formalism is presented, with 
particular emphasis on deriving the linearized field equations. In Section ^, we 
present a detailed discussion on the calculations of the shear viscosity for such higher- 
derivative theories. A thorough inspection of the field equations confirms that these 
existing methods are equivalent to reading off the shear viscosity from the coefficient 
of the kinetic term. Using this observation, we are then able to demonstrate explicitly 
that ghosts have infiltrated the earlier attempts at evaluating rj. Section ^ begins 
with a proposal for new EC's, as part of a revised prescription for evaluating the shear 
viscosity. This choice forces the decoupling of the ghosts that are inherent to higher- 
derivative theories. Subsequently, we explain how our new prescription protects the 
KSS bound and elaborate on the physical basis for the bound's validity. In Section ^ 
we resolve an apparent contradiction between our findings and the special case of 
Gauss-Bonnet gravity. Section ^ contains a brief summary of the results and their 
significance. The paper concludes with an appendix, where we provide a detailed 
description of previous methods for calculating rj. 

2. Equations of motion 

2.1 Specific 4-derivative models 

For clarity and concreteness, we consider the case of a 5-dimensional AdS black brane 
and focus on the class of four-derivative theories of pure gravity. Let us now set the 
notation and conventions. The background metric can be expressed as 

J 2 2 

ds^ = -F{r)dt^ + -^ + ^\dx^ + dy^ + dz^] , (2.1) 

r [r) L^ 



where F{r) = j^ 



1 - 1^ 



7^ 



while L and rh respectively denote the radius of curva- 
ture and black brane horizon. The function F vanishes on the horizon, F{rh) = . 
The theories of interest are completely described by the Lagrangian of Einstein 
gravity in AdS space, 

CE = n + ^, (2.2) 
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and three perturbative corrections as follows: ^ 

Ca = aL'^n^ , (2.3) 

Cb = PL'^TZa^U"-'' , (2.4) 

Cc = ^L^n^^cdn^"^' . (2.5) 

The validity of the gauge-gravity duality mandates the hierarchy rh ^ L ^ Ip 
(where Ip is the Planck length), thus making the dimensionless coefficients a, /3, 
7 ~ ip/L"^ much smaller than unity. We also consider, with particular emphasis 
in Section ^, the special Gauss-Bonnet combination of the three perturbations or 
Cgb = X[Ca- 4:Cb + Cc] . 

Our interest is in the case of weak gravity, and so the linear expansion of the 
metric about its background gab = Qab + ^ab + Olh"^] is valid. Similarly, we need 
linearized expressions for the various curvature tensors; for instance, 

T^abcd{h) = - [VcVfe/lad + ^d^ahbc " ^ d^ bhac " ^ c^ ahbd\ , (2.6) 

from which the contracted forms follow. An overline indicates the background ge- 
ometry. 

To obtain the linearized graviton field equations, we use the following identity 
0,131: 



- 2V„Vb 



5C 



}"^a b. 



+ n 



abcp 



5C 



S-^abc' 



- Igp^c = . (2.7) 



6gP'!i 

One can obtain the desired equations by expanding out each term to linear order 
in the /I's. (In our case, the first term is trivially vanishing.) For the current anal- 
ysis, this treatment is preferred over the more familiar procedure of varying with 
respect to the metric, as the absence of explicit derivatives in the Lagrangian means 
that boundary terms need not be stipulated. To determine t], we compute the xy 
component of Eq. ( p.7|) , knowing that the {x,y} sector does not mix with other po- 
larizations. The calculations are carried out in the transverse, traceless gauge and. 



^For a four-derivative theory, any derivative of the Riemann tensor reduces to a surface term 
and can therefore be ignored. 
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often, at the horizon. The brane horizon is the most appropriate surface for an anal- 
ysis of 1] from the gravity side; however, none of our conclusions are sensitive to this 
choice. Deriving the Einstein equation for Ce and multiplying through by a factor 
of —2, one then obtains 



Oe hxy = 



n^ + n + j^ 



I'xy 







(2. 



where Ul = U - Tl ^ -iVy - 2n\y and D = T^VaVb ■ The operator D^ is 
the spin-2 Lichnerowicz d'Alembertian [^ {i.e., the analogue of D for a graviton in 
curved space) for the background metric. When the field equation is re-expressed in 
terms of = g^^hxy, it is formally equivalent to the Klein-Gordon equation for a 
massless scalar throughout the spacetime. On the horizon, in particular, Eq. (|2.8| ) 
becomes 



n^ hxy — □ hxy — . 



(2.9) 



This observation proves to be important later on. 

Let us next consider the leading-order contributions from the corrections. After 
tedious but straightforward calculations, we obtain 

(2.10) 
(2.11) 

(2.12) 



Ca -^ -40q 


nhxy + --- , 






Cb ^ f3L' 


D^ - 2^^^nl 

r 


'^xy 1 ' ' ' ) 




Cc ^ 47^2 


V 2 


2 r r"^ J 


D 



'T'xy 



The ellipsis represent perturbative-order mass terms, which will be omitted for the 
remainder of the paper, as these are inconsequential to the current discussion. In 
any event, one can eliminate these by expressing the on-shell form of the full field 
equation in terms of = h^ . 

The n^ terms in the equation of motion are focal to the ensuing discussion, so 
let us be more explicit how they arise. After linearizing Eq. ( p.7|) , one finds that a 
D^ term can only come from V a^ b'^ {h)"" xy^ 5 where X"-'"^ = 6C/6Tlabcd ■ Con- 
sidering, for instance, the 4-index Riemann-squared case, one then has '^"'xy'' ~ 



2^L?'R^^y , which can be linearized via Eq. 



. In this case, only the second 
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term in Eq. ( p.6|) is relevant, leading to 27L V \nW ahxy ■ Next, commutator re- 
lations such as [Va,Vt,]V^ = T^abc ^d cau be utilized to attain (up to mass terms) 



27L2 



bac-z 



D^ + R^ VfeVc hxy or 27/.^ [D^ — -p-D] h^y , with the latter form follow- 
ing from symmetries of the background as explained below. Finally, restoring all 
numerical factors, one recovers the D^ term in Eq. (|2.12| ). 

An important caveat about Cc is that, for this case only, one finds terms that 
cannot be expressed directly in terms of D; with these always being of the form 
TZ r^V aSbhxy (as well as the obvious x ^ y analogue). However, we are interested 
in the t and r sectors because time derivatives are responsible for the ghosts (if any) 
and only radial derivatives will survive once the hydrodynamic limit is imposed to 



■tx 



tx 



calculate rj. Then, since only the terms with a,b = {r, t} are relevant and 7?. 
"^ rx = ~1^V' ' ^^ '^^^ make the substitution 7^ x'^a^bhxy — ^ —^^^Ohxy ■ 
The remaining discrepancy [^^""^^ ~ {~l^)] ^""^zhxy = \^[2F - rdrF] k'^hxy 
(with k = —iVz ) may be regarded as a perturbative-order mass term. 
The Gauss-Bonnet combination leads to the correction 

Cgb -^ -2XL^^n hxy , (2.13) 

where the background relations ^^ + 2^ = -^ and d'^F = 4^ — ^^ have been 
used to re-express the contributions from Ca and Cc respectively. As a useful check, 
we find that this Gauss-Bonnet correction reduces in the AdS vacuum {vh = 0) to 
—AXnhxy , in agreement with the known result PD| . 

For future reference, the horizon forms of Eqs. ( |2.1CI| - prT3| ) are as follows: 

(2.14) 
(2.15) 

(2.16) 
(2.17) 

The background metric has been used in the above and any subsequent compu- 
tations. Since our main interest is in corrections to r] up to leading order in ti/L'^-, 



Ca 


-^ -40an hxy , 




Cb 


^ (5L^ 


□= - ^□_ 


h-xy ; 


Cc 


-^ 4'yL^n^ hxy , 


CcB 


-^ -8A[ 


_l llxy ■ 





this choice is justified for contributions that come directly from the perturbative part 
of the Lagrangian. But what about those from the Einstein part? To understand 
why the background metric still suffices, let us consider the following: Given the Ein- 
stein Lagrangian and some perturbative correction, the leading-order effect on the 
near-horizon geometry can only be to shift L and Th from their background values. 
We can also, to leading order, re-express any contribution from the perturbative part 
in terms of the same shifted values of L and r^. Now, if our primary interest is to 
compute 7]/ s for a given theory, then all dependence on L and r^ will cancel out of 
this ratio. Meanwhile, if our main interest is to compare rj for two different theories, 
then it is natural to do so at a fixed value of brane temperature T = r/j/vrL^ and, 
hence, at common values of L and rh. So that, for current purposes, any correction 
to the background geometry never does come into play. 

We can restate the above argument in another way. As already discussed, one 
can fully determine r] for two-derivative gravity from the coefficient of the kinetic 
term in the Lagrangian or, equivalently, from the coefficient of the D terms in the 
field equation. For the Einstein part of the Lagrangian, the latter coefficient is 
universally —1/2. This is true anywhere in the spacetime and irrespective of the 
solution. So that, to understand how 77 is corrected for an extended theory, we 
need only to determine how the additional ( "non-Einstein" ) parts of the Lagrangian 
explicitly contribute. And, with these additions already being of order ti/L'^-, the 
Einstein background metric ( |2.1| ) suffices. 

2.2 Generic 4-derivative gravity 

Let us emphasize an important property of any four-derivative theory of pure gravity 
and, in fact, any higher-derivative extension that depends only on the Riemann 
tensor (and its derivatives): After the common factor of L? {cf, Eqs. ( p.3[ - p.5|) ) has 
been factored out of the corrected part of the field equation, the coefficient of the 
n^ term is a geometry-independent number (possibly zero) and, hence, a spacetime 
invariant. By geometry independent, we mean that it can depend on the form of the 
Lagrangian but not on the solution to the field equations. This outcome follows from 
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the observation that the prefactor for the D^ term is dimensionless; meaning that 
it can not depend on the Riemann tensor nor its derivatives, as any of these have a 
strictly positive mass dimension (nor can it depend exphcitly on the metric, as the 
Lagrangian is only a function of the curvature). On the other hand, the prefactor 
for the linear-D term has a dimension of mass squared, and so it can depend on the 
Riemann tensor. 

Consequently, for a generic four- derivative theory at any radius r, the corrected 
field equation can be written as 



^Ehxy + eL 






h,y = 0, (2.18) 



where e ~ Ip/L^ is a perturbative coefficient, a is a fixed number of order unity and 
b{r) is a radial function of order unity; all of which are dimensionless. On a constant- 
radius surface, b can also be regarded as a number. We are interested in corrections 
to r]/s to leading order in e and therefore need to consider a finite (small) e. 

We need not limit considerations to the {x, yj-polarization channel for the gravi- 
tons, as 7] can also be extracted from the so-called shear and sound channels (see, 
e.g., Appendix B of [0). For the shear channel (the sound channel follows along 
similar lines), the relevant gravitons are hra, hta, hza with a = {x,y} . The field 
equations for this class are more complicated but can be simplified by choosing the 
radial gauge hr^ = { fi = {r,t,z,a} ) and then considering the gauge-invariant 
combination Z = g^^khtx + g^^oohzx ( uj = +i'Vt , k = —iVz )• In this way, the 
coupled set of equations is turned into a single equation for a single graviton mode 
Z. This equation is necessarily of the same generic form as Eq. ( p.l8| ), as none of 
the above arguments are specific to the {xyj-polarization class. Then, just as we go 
on to show how the shear viscosity can be extracted from the kinetic coefficient of 
hxy, the same can be deduced for the kinetic coefficient of Z |^ (and analogously 
for the sound channel 
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3. Calculating 77 

3.1 Contributions to r] 

Next, we consider the shear viscosity rjx for the extended theories Ce+jC-x , with X = 
{A,B,C,GB}. But, first, let us recall what is known about two-derivative theories: 
The kinetic coefficient of the transverse gravitons at the horizon fully determines the 
shear viscosity for both the "graviton fluid" on the black brane and its field-theory 
dual at the AdS boundary. This follows from a rigorous calculation in ||31| (also, ||43|| ) 
that applies to any two-derivative (but otherwise arbitrary) 5D gravity theory. There 
it was shown that the Kubo form of the field-theory viscosity tjft can be expressed 
purely in terms of the bulk geometry. They then imposed the standard choice of 
incoming plane wave BC at the horizon and the Dirichlet BC at the outer boundary 
to arrive at 

1 T 
VFT = Vgrav = J^^^i^ = ^h) , (3.1) 

where K is the kinetic coefficient of the transverse gravitons (here, normalized so 
that its Einstein value is unity). 

Supposedly, we can determine how each extension modifies the Einstein viscosity 
t]e by inspecting the horizon forms of the corrected field equations ( |2.14[ - PTT7D . The 
kinetic coefficients, in particular, should tell us how the viscosity is modified. The 
scalar-squared and Gauss-Bonnet cases {A and GB) are effectively two- derivative 
theories, and so the associated corrections can be read off directly from Eqs. ( ^.14] ) 
and ( p.l7|) . We then have rjA = tjeII — 40a;] and tjgb = Ve[^ — 8A] ; both in 



agreement with the already-known results. ^ 

However, for the Ricci and Riemann-tensor squared cases [B and C), there are 
now n^ terms present in the field equations. Following the standard prescription, one 
would conclude from Eqs. ( |2.15| ) and ( |2.1(j| ) that ris = ?7£; [1 — 8/3] and f]c = Ve , 



^For rjGB, see (e.g.) M. That the scalar-squared case (A) is a match, even if not documented 
exphcitly, follows from the confirmation of tja/sa — Ve/se , as must be true for any f{TZ) theory. 
Note that sa = sb [1 - 40a] follows from Wald's formula ^ ^. 
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which are indeed in perfect agreement with those obtained from any of the previously 
described methods. ^ 

On the horizon, dh^y and the background field equation OEhxy are inter- 
changeable up to order Ip/L'^- Then, as any D must be of order Ip/L"^, D^ ^ D 
certainly follows. Off the horizon the situation changes, since it is \I\Ehxy = 
[n + 2F{r)]hxy that vanishes on shell, and so D^ ^ D is no longer valid. A 
single n^ term now makes a generic contribution of 

D^ = n| - AFDe + ■■■ , (3.2) 

where the dots signify that we are disregarding mass terms as usual. In terms of our 



generic theory of Eq. ( 2.18 ), 6 — )■ 6 — AaF . The D^ term has, through its prefactor 
a, entered into the calculation of the kinetic coefficient. ^ 

It would seem that the horizon value of b and, hence, t] are insensitive to the 
n^ terms. We will, contrary to this expectation, show that the D^ terms are always 
implicated with ghosts; with these directly influencing the value of b throughout the 
spacetime — including at the horizon! 

3.2 Ghosts in the machinery 

It is well known that the presence of a D^ in the equation of motion signals that a 
ghost lurks in the theory. The usual argument for ignoring this is that the Planck- 
scale ghost becomes infinitely massive once the ultraviolet cutoff has been sent to 
infinity [H|. However, for the case at hand, e ~ ^p/L'^ needs to be kept finite. If 



one takes the strict decoupling limit e = and all corrections originating from the 
higher- derivative terms vanish, as expected. 



^Similar to case A, the Ricci-square result follows from knowledge of the entropy density via 
Wald's formula and that tje/sb = tje/se R1- Also, to avoid any confusion, the ratio rjc / sc 
changes by an overall factor of 1 — 87 + ©[7^] , same as for Gauss-Bonnet. Only that, for the 



Riemann-squared case, this correction comes entirely from the entropy |29 . 

^This distinction between the horizon and other surfaces seems contradictory to the claim that 

rjFT can be calculated at any radius in the spacetime |30[| . However, this claim has only been 

established for two-derivative theories. 
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Let us begin here by exposing the ghosts. We will do so using a generic four- 
derivative theory, as described by Eq. ( p.l8| ), to emphasize that the ghosts are neither 



model nor radius specific. First, Eq. ( |2.18|) can be recast, using Eq. (|3l^ ), into 

Ue + eaL^Dl + efoD^ h^y = , (3.3) 

where b = b — 4aF is a radially dependent parameter that generally differs from b 
but a is the exact same number as before. Or, up to leading order in e, 

(1 + eb)DE [I + eaL'^DE] Ky = , (3.4) 

with the tilde on b henceforth implied. 

The factorization in Eq. ( |3.4| ) implies that, at order e, there are two separable 
modes in this problem. Let us then write the graviton as H = Hq + Hi (with tensor 
indices implied), where these are meant (and shortly will be shown) to satisfy 

□e^o = , (3.5) 

1 



"'°+L2 



Hi = . (3.6) 

To explicitly expose the two distinct modes, let us invoke the standard practice 
in flat spacetime of inverting the linearized operator and then factorizing. That the 
inversion process extends in a straightforward manner from flat space to AdS space 
follows from the fact that D^; is a scalar. Then D = —p"^ of flat space is replaced 
with an appropriately generalized momentum; say D^; = —p ^ . On the horizon of 
a black brane, p^ = p^ . Then, 

1 = 1 1 (37) 

Ue[1 + eaL^UE] Ue Ue + (eaL^)-! ^ ' ^ 

The flrst term on the right can be identifled as the propagator of the massless Einstein 
graviton. Meanwhile, the second term is the propagator of a massive ghost graviton, 
as evident from the "wrong" sign relative to the former. This outcome does not 
depend on the sign nor magnitude of the perturbative coefficient e, which only de- 
termines the mass of the ghost. Let us recall that a,6 ~ (9[1] and e ~ I'i/L'^ , and 
so the ghosts have a squared mass M^ ~ {eL"^)^^ ~ 1//^. As already stressed at the 
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end of Section 2, we are looking at corrections to leading order in e. Consequently, 
the very massive ghosts may indeed affect the result. 

If a four- derivative theory does have ghosts, these persist throughout the entire 
spacetime because of the invariance of the dimensionless constant a associated with 
the n^ term. ^ Hence, the ghosts are not confined to the AdS bulk interior and can 
contribute to the field-theory viscosity to leading order in e. 

Now suppose that we want to restore the normalization factor of (1 + 6) to the 
mode equations ( p.5| , |3.6| ). According to any of the standard methods, the factor 
(1 + e6) is strictly attributed to the massless graviton. This follows from rj always 
going as the value of (1 + eh) at the horizon (see Subsection 3.1) and the implicit 
assumption that any ghost mode would have ultimately decoupled. On this basis, 
the mode equations should be reformulated as 



taU + — 



(1 + eh)UEn^ = , (3.8) 

Hi = . (3.9) 

But, as 6 = b{r) does not commute with D^;, this formulation is problematic in that 

it could not be obtained through a process of factorization (c/, Eq. ( p^) ). Hence, 

whatever might be the actual coupling for the massless graviton, it can be generically 

different than the standard result of (1 + e6). This difference can be attributed to 

the fact that the gravitational coupling (l + eb) can only be attributed to the sum of 

the modes and not to any single one of them. This argument shows that the ghost 

mode can indeed impact the calculation of t]. 

In summary, we have seen that the ghosts have, through their association with 

the coupling correction b, managed to infiltrate into the coefficient of the kinetic 

terms and, thus, into the calculation of r]. This is true insofar as the squared mass 

in units of the AdS radius L^M^ ~ e^^ ~ -^^/^p is kept large but finite. In the next 

section, we will propose a procedure that guarantees the complete decoupling of the 

ghost mode and allows for an unambiguous calculation of rj/s. 

^The same claim can, in fact, be made for any higher-derivative theory that depends only the 
Riemann tensor. The inclusion of other types of fields could jeopardize the invariance of a, but this 
parameter would still be typically non- vanishing at all radii. 
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4. New (and improved) prescription for calculating rj 

We conclude that, to proceed, one should ensure that the influence of the ghosts 
is completely eliminated before the calculation of t] is carried out. This could be 
accomplished through a carefully imposed choice of BC's, as we explain next. 

4.1 An exorcizing boundary condition 

Let us reconsider the generic field equation ( p.4|) . For a generic extension of Einstein's 
theory, one can expect two BC's per each order of e. This enables one to individually 
fix each order of the solution to fulfill the stated conditions [^. Here, we have 
rearranged the solution so that it is separated according to the different degrees of 
freedom rather than perturbative order. Nevertheless, as is evident from Eqs. ( |3.5| ) 
and ( |3.6| ), both of the modes (massless graviton and ghost) still satisfy a simple 
quadratic equation and, thus, both have a corresponding incoming and outgoing 
solution. And so, with four BC's at our disposal, we are using one of these to 
stipulate that the ghost is an incoming plane wave at the horizon and another to 
kill off the ghost by, for instance, imposing the Dirichlet BC at the horizon. ^ The 
remaining two BC's are then used to ensure that the massless graviton satisfies the 
usual pair of BC's: Dirichlet (or vanishing wavefunction) at the AdS boundary and 
incoming plane wave at the black brane horizon. 

Our new BC's are compatible with the standard set of incoming plane wave at 
the horizon and Dirichlet on a radial shell ^ for the total wavefunction. Let us recall 
that the total solution can be expressed as H = Hq + Hi , with Hq/i respectively 
labeling the massless graviton and massive ghost. We also recall the relevant pair of 
quadratic equations, 

DeHo = , (4.1) 



^Generically, these hydrodynamic modes are regular everywhere, except at the horizon where 

regularity must be imposed by hand (see, e.g., pq]). Hence, fixing the wavefunction to vanish at 

the horizon is a sufficient condition for the mode to be vanquished throughout the spacetime. 
®As stressed in [M, there is nothing special about the outer AdS boundary in this regard, even 

if the standard choice. The Dirichlet condition can readily be imposed on any radial shell exterior 

to the horizon. This includes at the horizon when taken as a suitable limit of the stretched horizon. 
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eaDE + -^ 



H, =0 , (4.2) 



with both admitting plane-wave solutions. 

Adopting the standard incoming plane- wave solution and using the dimensionless 
(inverted) radial coordinate u = rf^/r"^ , we can express the total wavefunction of 
the mixed-index gravitons $0/1 = Q^^Hqji as the sum $(u,a;,g) = Co$o + Ci<l>i 
such that 

•^0 = [f{u)r-"To{u,uj,q), (4.3) 

$1 = [/(«)]-^"/'Ti(«,a;,g). (4.4) 

Here, f{u) = l — v? (essentially, \gtt\), ^ and q represent a dimensionless frequency 
and wavenumber (as defined in, e.g., [0), Cq/i are normalization constants and Tq/i 
are mo del- dependent functions of u,u,q . The latter functions can be uniquely fixed, 
up to normalization, by imposing regularity at the horizon (with regularity assured 
elsewhere). 

The hydrodynamic limit implies a;,g ^ 1 , so that — in practice — one expands 
$0/1 out in increasing powers of u and q and considers only the first few terms in the 
series. For instance, Einstein gravity ( e = ) leads to Tq = 1 + 0[uj'^,q^] ^^. ^° 
Technically, however, the hydrodynamic limit should only be applied at the very end 
of the calculation. 

For $0 ; the Dirichlet BC is imposed at the outer boundary ( m — t- ) and then 
the dispersion relation u = uj{q) is fixed to attain 

[/(n)]— /2To(«,u;,g)L^„ = 0. (4.5) 

Meanwhile, the normalization Cq is fixed in the standard way 

C,' = [/(n)]-^-/'To(^,a;,g)|^^^ , (4.6) 



^'^At the horizon, either of our modes has precisely this same form. For the massless graviton, 
this is trivially so at any radius. For the massive ghost, this follows from the red shift at the horizon 
making this mode effectively massless. 
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where m is a constant "cutoff" scale. Consequently, the renormalized wavefunction 
$0 = Co$o satisfies 



$r 



1 



u-i>0 



(4.7) 



For the ghost mode $i, we now deviate from the usual convention and rather 
fix its normalization with the constraint Ci = 0. Since regularity has already been 
imposed at the horizon through the correct choice of the functional form Ti, the 
renormalized wavefunction $i = Ci$i is compliant with the horizon Dirichlet BC 



Hi 







u— )-l 



(4. 



As a result, the total normalized wavefuntion $ = $o + *^i reduces to $0; 



^{u,u,q) 






-1 -iu)/2 



To{u,uj,q) 



To{u,uj,q) 



(4.9) 



M^O 



In this way, we have ensured that the amplitude of the heavy ghost mode is made 
to vanish without insisting on the vanishing of e. 

The effect of imposing the BC's that "exorcise" the ghost is to set the right-most 
term in Eq. ( |3.7] ) to vanish. We can then re-invert the remaining operator to obtain 
a "unitarized" version of the field equation ( p. 41 ) in which the rightmost factor has 
been set to unity: 

{l + ec)nEKy = 0, (4.10) 

where the order-unity radial function c{r) is related to but not necessarily equal to 
the original radial function b{r). As already discussed, this normalization factor may 
well have changed in the described process. 

Because the value of the function c could have been "infected" by the ghost 
mode, it is desirable to bypass the individual evaluation of r] and s, which would 
require an exact calculation of c. Rather, we will calculate the ratio r]/s from which 
the exact value of the normalization function c actually drops. 

The field equation Eq. ( [4.1CI| ) can be derived from an effective action with the 
following Lagrangian: 



Ceff = - 2 (1 + ec) r'V.h^y y,Ky + 



(4.11) 
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Since this is now a two-derivative theory of gravity, we can directly invoke Eq. ( 3.1 
to determine the shear viscosity; hence, given c, rj goes as (1 + ec)rjE ■ 

Now, what about the entropy density? Once the ghosts have been exorcized, 
the gravitational theory is Einstein's with a renormalized gravitational coupling. 
This follows from |^^ ^, ^, where it has been made clear that a consistent two- 
derivative theory for a massless spin-two field is uniquely described by Einstein grav- 
ity. In general, differently polarized gravitons would "perceive" different values for 
the couplings. However, for a consistent two- derivative theory, there can be no such 



distinction between the polarizations |^, ^, ^. So that, by eliminating the ghosts, 
we have also rendered gravity to be polarization indifferent. 

We can then call upon the analysis from Section III of [^, where the Wald 
entropy has been reinterpreted as an effective coupling for the r, t-polarized gravitons, 
to deduce that s goes as (1 + ec)sE ■ Meaning that the ratio rj/s must necessarily 
be equal to the Einstein value of 1/Att. That is, the "ghost-busted" theory may have 
a shifted value of r], but it comes with the assurance that s will be corrected from 
its Einstein value by the very same amount. 

4.2 Comparison to another choice of boundary conditions 

We would like to be precise on how our newly proposed prescription differs from those 



previously used; in particular, the effective-action method ||33| and those similar to it 
{e.g.) [^. Given a four- derivative theory, a prescription for calculating r^ is a three- 
step process: (i) find the equations of motion, (ii) solve the equations of motion, (iii) 
decide on how to extract the viscosity. There obviously can be no dispute on the 
first step and, given that the action has been consistently reduced to a two-derivative 
form, there should be no contention about the third. Where our new prescription 
then differs is only in the second step, where a clear distinction has been made 
between the appropriate choice of EC's and, hence, on the resulting solution. It is 
interesting to compare these two sets of solutions in detail, as we do next. 



Let us begin here by recalling Eqs. ( |4.10| ) and ( |4.2| ); these being the field equa- 



tions for the (post-normalized) massless graviton Hq and the ghost graviton Hi, 
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respectively: 

{1 + c)UeHo = 0, (4.12) 

1 



eaDE +2 



Hi =0 . (4.13) 



We want to compare these equations with those obtained in the "usual" manner 
by which one expands the solution according to perturbative order (rather than 
isolating the degrees of freedom). Let us label the zeroth-order and first-order modes 
respectively as Hq and ehi. Then, the previous field equations translate into 

OeHo = , (4.14) 

[eaL^DE + eb] DeHo + eD^/ii = . (4.15) 



The standard choice of BC's |^ is incoming at the horizon and Dirichlet at the 
outer boundary for both Hq and hi. 

To avoid clutter, let us consider the simple but instructive case of b = c . 
Comparing the two sets of field equations, one finds that compatibility (to first- 
perturbative order) requires the following dictionary: 

ho = Ho + ^Hi, (4.16) 

hi = -aL^DEHo - ^bHi , (4.17) 

L-Ib 

so that the two sets of modes share a non-local relationship. Now, it follows from 
the above equations and D^; ~ e (so ^ ~ ^ ) that Hq is ho with a perturbative 
correction whereas Hi is strictly is of order e, as it must be since the ghost is required 
to vanish in the Einstein or e — )■ limit. 

Comparing our exorcizing BC's to the standard set, we see that the difference 
between them is quite subtle, as it only appears at order e for both /iq and hi. From 
the discussion in the preceding subsections, it is clear that the values of rj and of rj/s 
are sensitive to the choice of BC's. In particular, an order-e difference in the BC's 
for the massless graviton results in an equally small but still significant change in 
the value of rj/s. 
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4.3 General arguments 

For theories leading to equations of motion with more than four derivatives, it is 
quite possible that not all of the additional degrees of freedom are ghosts (see the 
discussion on six- derivative equations at the end of this section); meaning that the 
ghost-reduced action need not be limited to two derivatives nor even be local. The 
latter because our prescription would, in general, yield a non-local field equation due 
to summing over two or more inverted operators. 

For such cases, we propose proceeding as follows: First, after re-inverting the sum 
of the non-ghost propagators, one should convert this ghost-free field equation into a 
local form by expanding the propagators of massive modes (with masses Mi) in terms 
of O/Mf ^ 1 . Then, following the effective-action method I^S], one should reduce 
the order of this field equation by treating the higher-derivative terms perturbatively 
while iteratively applying a lower-order field equation and the BC's. However, the 
initial field equation for starting the iterative purposes is, in general, no longer that of 
the background. It is, rather, the equation one obtains by perturbatively expanding 
and then suitably truncating the ghost-free field equation. Finally, one should, as 
before, translate the reduced field equation into a two-derivative effective action and 
then apply Eq. (|3.lD. 

One might then wonder as to how rj/s would turn out for such higher-derivative 
theories. We will argue below that, on general grounds, this ratio can now change 
but only in such a way that it increases relative to its Einstein value 



We have advocated elsewhere |T^ that any gravitational coupling, such as r], 
can only increase from its Einstein value for a unitary extended theory. This is 
because any consistent extension will necessarily introduce new degrees of freedom 
to supplement the Einstein graviton. Given that these are unitary, the extended 
theory can only introduce additional channels that act to increase the couplings. 

Let us recall the relevant discussion from [^. Given our gauge choice, the 
hxy gravitons can only couple linearly to other particles of spin 2 and of the same 
polarization. Near the horizon, the associated IPI graviton propagator ^^ takes on a 



11 



See (e.^.) [18J for the fully decomposed form of the graviton propagator. Note that this flat- 
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particularly simple form: 

i ^ 

where q^ = —q^q^i is the spacelike momentum, pE is the gravitational coupling 
for the Einstein graviton, and the p^'s are the couplings for the additional spin-2 
particles of an extended theory — any of which could be massive ( mf 7^ ) or 
massless ( m^ = ). ^^ The couplings can, as indicated, depend on the energy scale 
q ; and we work in units such that P£;(0) = 1 correctly fixes the Newtonian force at 
large distances. 

The essential point here is that the couplings are really spectral densities of the 
schematic form p = ^„(0|/i|n)(n|/i|0) [^ and, as such, assured to be positive as 



long as all the inserted states have a positive norm. Meanwhile, the Einstein coupling 
can itself be modified in either direction (c/, Eq. ( [4.10D ), but this is an illusionary 



effect, as one would always recalibrate her instruments to maintain the p_b(0) = 1 
normalization. In any event, any such modification to the Einstein coupling would 
immediately cancel out of the ratio rj/s, as Einstein gravity must be insensitive to 
the polarization of gravitons (which is just a restating of the equivalence principle). 
Conversely, the "non-Einstein" gravitons are generally sensitive to the polarization 
but can still only act to increase t]/s from the Einstein value of l/An. 

To sum up, the microscopic theory is, on the basis of unitarity, telling us that 
T]/s can only increase from the Einstein value. This point allows us to clarify what 
was previously meant by "considerations of unitarity" . For addressing such matters 
of principle, a Planck-scale ghost is obviously unreliable, as would be the case for 
any mode near or above the effective-theory cutoff. On the other hand, including 



space form of the propagator is appropriate near any horizon — including that of a brane in AdS 

space — as any horizon has an effectively flat geometry. This is so because D^ (x —df + d^^ as 

r — > r/i (where r» = j dr/F is a generalized "Tortoise" coordinate). 

^^By the same reasoning as in the prior footnote, any particle near the horizon is effectively 

massless and "perceives" an effective 2D geometry. Hence, one can set all the mf to zero and need 

only consider the t and r» components of q^. Note, though, that the full 5D tensorial structure of 

the gravitons is still maintained. See |50[| for further explanation. 
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the ghost in the calculation enables r]/ s to change freely in any direction. Hence, 
we conclude that, to obtain results from the effective theory that are consistent with 
microscopic unitarity, it is necessary to impose EC's that eliminate these unreliable 
modes. Consistency then requires that a physically meaningful quantity such as t]/ s 
should no longer be sensitive to cutoff-scale physics. 

Applying the preceding ideas to our generic four-derivative model, we expect 
that the ratio r]/ s must remain independent of the perturbative coefficient e at linear 
order. If this were not so, one could always reverse the direction of the t]/ s correction 
by simply changing the sign of e. Reassuringly, this is exactly what our revised 
prescription ensures! After all, the non-Einstein degree of freedom is inevitably a 
ghost (c/, Eq. ( p.Tp ); so that, with our choice of EC's, the physically relevant theory 
is guaranteed to reduce to Einstein's. 



4.4 Six and higher-derivative equations of motion 

Let us now discuss how the situation that we have described in such detail for the 
case of four-derivative theories is modified for theories leading to six-derivative equa- 
tions of motion. We will only provide a sketch, relegating a detailed discussion to 
a subsequent paper. Six or higher- derivative equations can be quite different from 
four- derivative theories. With six derivatives, for example, one would obtain two 
additional degrees of freedom, only one of which need be a ghost. To illustrate this, 
let us consider "toy-model" versions of a six-derivative linearized operator. For in- 
stance, Qi = De + Cni (c ~ -^j ~ C^n|; (with ( a perturbative coefficient, c y^ 1 
parameterizing mass and units of L = 1), which factorizes to give 



(4.19) 



In this case, both of the massive gravitons are ghosts and should be eradicated from 
the calculation. On the other hand, Q2 = De — (^"e l*^ + c) "'" C^^^e leads to 

^2 
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(4.20) 
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This second case is more interesting because only one of the massive gravitons 
is a ghost, with the choice depending on whether c is greater or less than unity. For 
such a model, it can then be expected that r]/ s does change, as a new unitary degree 
of freedom would generally be sensitive to graviton polarization. Nonetheless, seeing 
that Q2 only differs subtly from Qi, one is now able to envision a scenario whereby the 
emergence of a new degree of freedom is correlated with a strictly positive correction 
to 7]/ s. Provided that the ghosts have been properly handled, it is our contention 
that this is, indeed, what must happen. ^^ 

5. Gauss— Bonnet redux 

There is still an important example to be dealt with; the Gauss-Bonnet gravity 
model Ce + Ccb ■ This theory leads to a second, rather than fourth, order equation 
for the graviton as can be seen from Eq. ( p. 171) . Hence, the issue of the existence 



of ghosts or choice of boundary conditions becomes irrelevant. Obviously, unitarity 
considerations cannot be applied in a useful way in this case. 
Yet, we have found (as did past studies such as [0) that 

^ = -^[l-8A], (5.1) 

Sgb 47r 

as Sgb is known to be functionally equivalent to se- Evidently, the ratio will decrease 
below l/47r whenever A is positive. So, this theory seems to provide a counter- 
example to the main theme of our paper, as it apparently demonstrates a ghost-free 
theory for which the viscosity-entropy ratio can be less than its Einstein value. 
However, this conflict is only an apparent one, as we now explain. 

There are two distinct ways of viewing the Lagrangian Ce+Cgb ; as representing 
a theory unto itself {i.e., distinct from Einstein's) or as a perturbatively corrected 
version of the Einstein Lagrangian. If one chooses the former point of view, then the 



"'^■^Onc might be concerned that s — also being a type of gravitational coupling — should likewise 
increase, possibly faster than 77. Nonetheless, we have shown elsewhere |15[ that it is always possible 
to choose fields and coordinates such that s is calibrated to its Einstein value. Then, since rj > r]E 
is a "gauge" -invariant statement, the previous claim follows. 
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microscopic theory is Gauss-Bonnet gravity and there can be no contradiction with 
what we have aheady said. All methods of calculations agree as they must, and the 
issue of imposing boundary conditions in a higher-derivative theory never arises. 

On the other hand, from the effective field-theory perspective, Einstein's theory is 
viewed as the infrared limit of the quantum gravity theory and the various corrections 
(£x) are regarded as a consequence of integrating out heavy degrees of freedom. 
Then, even if the Lagrangian has the finely tuned Gauss-Bonnet form ( /3 = —4a = 
—47 ) at some scale, it will be of the generic form at other scales and lead to quartic 
equations of motion for the graviton. 

So the issue now becomes an order of limits: Should one first impose the Gauss- 
Bonnet fine-tuning condition and then calculate r^/s, leading to Eq. ( ^.1|) and possibly 
to values of 7]/ s lower than l/47r? Or should one first eliminate the ghosts from Cb-, 
Lci then impose the Gauss-Bonnet fine-tuning condition and only afterwards evalu- 
ate the ratio? If one has a reason to regard the Gauss-Bonnet theory as fundamental, 
then the former choice is correct. With this interpretation, Gauss-Bonnet gravity 
can no longer be regarded as a perturbative correction of Einstein's theory and, so, is 
outside of the scope of our previous claims. On the other hand, if the Gauss-Bonnet 
model is to be regarded as the leading-order modification of the Einstein Lagrangian, 
the latter choice should be made and the correct calculation leads to a renormalized 
Einstein theory with the outcome tjgb/ Sgb = l/47r . 

6. Concluding discussion 

Let us summarize: We have shown that, for higher-derivative theories that extend 
Einstein gravity, earlier choices of BC's for the purpose of calculating the shear vis- 
cosity have resulted in contributions from ghosts to the value of 77. We have proposed 
a different choice of BC's that decouples the ghosts. An immediate consequence of 
our proposal is that, when the extensions are limited to four derivatives, the ratio of 
shear viscosity to entropy density saturates its Einstein value to leading order in the 
strength of the corrections; that is, saturates the KSS bound. We have gone on to 
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argue, on general grounds related to unitarity, that our prescription can be expected 
to enforce the KSS bound for theories leading to six and higher- derivative equations 
of motion and to all orders in the strength of the corrections. In such general cases, 
the saturation of the bound should no longer persist, as any gravitational coupling 
would naturally increase (from Einstein's) after additional unitary degrees of freedom 
have been introduced into the gravitational sector. The validity of our claims under 
these more generic circumstances remains an interesting, open question, which we 
hope to address with explicit calculations. 
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A. Previous proposals 

Let us be more specific on how the shear viscosity has previously been evaluated 



for a higher-derivative theory. We begin with the "effective-action" method |33] and 
remind the reader that other approaches — in particular, the generalized canonical- 
momentum treatment of [^ — are know to yield identical results. 

Let us consider the Lagrangian £ = Ce + ^Cx , where Cx is an arbitrary four- 
derivative extension and e controls the perturbation. To quadratic order in the xy 
gravitons of the metric, C can be written as 

^^C = A{u)(j)"(f) + B{u)(f)'(j)' + E{u)(f)"(f)" + F{u)(f)"(f)' + ■•• , (A.l) 
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where terms of lower- derivative order, as well as terms vanishing in the hydrodynamic 
limit of vanishing frequency and transverse momentum, are denoted by "■ ■ ■ ". In this 
appendix, we use the radial coordinate u = 1/r^ and set L = r^ = I so that the 
horizon is at Uh = 1 ■ Additionally, = h^^ , a prime denotes a differentiation with 
respect to u and A,B,E,F indicate mo del- dependent radial functions. (For Einstein 
gravity, E = F = .) The background metric ( p.l|) becomes 



ds' 



m,e 



du" 



H [dx"^ + dy"^ + dz^] 



(A.2) 



u Au'^f{u) u 

with f{u) = 1-^2 and /(I) = . 

As a first step, one integrates by parts (with the relevant surface terms presumed 
to exist \E^) to obtain 



^/^C = E(t)"<p" -f 

The associated field equation is then 

F'' 



B-A 



F' 



+ 



(A.3) 



E(t)" 



+ 



A-B 



+ 



A-B + 



F' 



b' + 



0. 



(A.4) 



Next, one is instructed to use the zeroth-order or Einstein field equation, cj)" = 
—Z'(j)' + ■ ■ ■ with Z = [ln{y/ —'gg^"')] , to reduce the higher-derivative E term, as 
E is already of order e. Twice applying the zeroth-order equation in the prescribed 
manner, one ends up with 

F'' 



(ZE)' - Z^E + B-A 







(A.5) 



and the omitted terms now also include those of order e^. 

Eq. (|A.5|) implies an effective Lagrangian of the following form: 



-gC 



eff 



r 



(ZE)' - Z^E + B-A- 



r^'du 



+ 



(A.6) 



The shear viscosity can be identified with the coefficient of the kinetic term for the 



Vs at the horizon |29|, |30|, pT|, and so one would deduce that 



n oc 



{ZE)' - Z^E + B-A- 



F' 



(A.7) 
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where the horizon hmit is imposed only at the end. 

Let us confirm that this result is indeed equivalent to that obtained from the 
canonical-momentum formalism of Myers ei. al. |^. Generalizing the methodology 



of [0 , these authors have essentially identified rj with the (inverse) ratio of to its 
canonical conjugate, as computed at the AdS boundary and in the hydrodynamic 
limit. A nice consequence of the latter limit is that the same identification can be 
made on any radial surface in the spacetime, including at the black brane horizon. In 
this manner, for a generic four- derivative theory, the shear viscosity has been formu- 
lated strictly in terms of the horizon values of A,B,E,F and the metric. Retrieving 



their formula from Eq. (3.22) of [Q, we obtain (after some minor manipulations) 
^ « -=L Y\i^y'\ +B-A-- , (A.^ 



where Y = \/—gu'g^^ ■ Comparing the last two equations, one can see that it is 



E ^^/] 



on 



enough to establish the equivalence of £i = {ZE)' — Z'^E and £2 = y 

the horizon. For this purpose, it is useful to recall that E is the coefficient of the 

four-derivative term (c/, Eq. (A.l)); meaning that it generically contains two factors 



of ^" ~ Y and is quadratically vanishing on the horizon. Hence, we can rather 
work with W = E/ p { f r^ Y ), with W assured to be regular as n — ?■ 1 . 

Let us now specialize to the background metric of Eq. ( [A.2| ), although, as one 
can check, this is not a necessary requirement. Making the substitutions E = pW , 
Z = y — i , Y = 2y/uf into £1 and £2-, one finds after simplifying that 

ifi = f'W + [f' - W , (A.9) 

Finally, after imposing the horizon limit, one ends up with 

j£,^2 ^ ni)w{i) + ni)w\i) (A.ii) 

in either case. 

It is instructive to focus on a specific gravity theory. Let us consider, for instance, 
the Riemann (tensor) squared model or C = Ce + Cc = TZ + 12 + •jTZabcd^"''^'^'^ , ior 
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which the background metric ( |A.2|) is appropriate and the exphcit forms of A,B,E,F 
are aheady given in Eq. (3.22) of |p6|. ^* Substituting these into Eq. ([A.7|) or 



Eq. (|A.8|) and then simphfying, we have 

V ^ 1 + 87 [/'(I) -/"(I)] , (A.12) 

where the Einstein term has been "normahzed" to unity and note that the correction 
term is, for this particular background, a vanishing quantity. 

Let us next consider the approach where one continues to identify rj with the 
horizon coefficient of the kinetic term for the xy-polarized gravitons in the original 
action. The premise being that the gravitational coupling for a given class of gravi- 
tons can be extracted from their kinetic terms |^ and that the shear viscosity is 



a measure of the coupling for the xy gravitons. One can, in principle, make this 
identification by expanding out the Lagrangian to quadratic order in h^y. It is, how- 
ever, simpler to extract the kinetic term by a method proposed in [^, which can be 
viewed as a generalization of Wald's entropy formula |^, ^ ^ and is tantamount 
to varying the Lagrangian by TZ^^^y . For the just-discussed Riemann-squared model, 
this implies 

r/ oc 1 + 47 \nZy]^^_, = 1 - 47/(1) . (A.13) 

Here, the correction term is also vanishing but the functional form apparently differs 
from that of the previous computations. Nevertheless, we can help to establish the 



claimed equivalence (for four-derivative theories) as follows: 

'V _ _ 

xy '^ X '^ xu '^ xt 



First, let us expand the curvature component out as 7^ =1Z x~^ xu~T^ xt~ 



IZ ^2 . Next, we use the equality IZ ^ = IZ f- , which follows from the spacetime 
being static and Poincare invariance on the brane. Making this substitution and 



■xy 



xy 



then expanding out IZ ^ in terms of its four-index constituents, we have TZ 
TZ tu'^ "^ ty ~^ T^ tz ~ T^ xu ~ '^ xz ■ Straightforward evaluation at the horizon 
then leads to IZ ' ^y = 2 [/'(I) — /"(I)] • Inserting this into Eq. ( [A.13| ), we verify a 
functional form for rj that agrees perfectly with Eq. (|A.12| ). 



-'^^Their ci is the same as our 7 and they denote A,B,E,F without tildes. Otherwise, the conven- 
tions are basically in agreement. 
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Such agreement does not generally persist for theories with six or more deriva- 
tives. For instance, let us consider the six-derivative corrected theory L = Ce + 



(T^abcdJ^ '^"'e ' fo^ which the effective-action method has already been used in [52 



to obtain rj = r]E[l — 32(^] . On the other hand, extracting the shear viscosity from 



the real action as prescribed in [^, one rather finds that rj = rjE [1 — 16^] • Note 
that the equations of motion for this model are limited to terms with no more than 
four derivatives, so that the two methods can already be in disagreement for theories 
leading to four-derivative equations of motion. 
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